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1. Introduction 

The original Gribov-Zwanziger confinement scenario [[T], ||] predicts an IR vanishing static 
gluon propagator Z)( |/?| ) in Coulomb gauge. The gluon propagator is also at the heart of varia- 
tional estimates to the ground state wave function [0, H]. A cross check of continuum results with 
lattice calculations had, however, long failed. A first study for SU(2) at fixed /3 = 2.2 indicated 
compatibility with Gribov's formula in the IR but was inconclusive in the UV [Q]. Later studies 
in SU(2) and SU(3) showed for strong scaling violations and a UV behaviour at odds with 

simple dimensional arguments |p. All these works calculate fixing the Coulomb gauge 

only at a given time-slice, neglecting the residual gauge freedom of temporal links. Also one takes 
for granted that multiplicative renormalizability holds for the full propagator D(\p\,po), although 
perturbative results point at a more complex picture [g]. We report here the results first obtained 
in JTo|], where a different strategy was adopted, fixing analytically the residual gauge and then 
studying the renormalization of the full spatial propagator D(\p\,po). We show that the Gribov 
formula [|l|] perfectly describes the lattice data for the static propagator. We refer to our original 
paper JTOj ] for details about conventions and notations. To achieve a good gauge fixing we adapt 



the algorithms developed in [11, 12]. 
2. Results 

The first observation made in [ |l0| ] is that the lattice bare propagator Dp(\p\,po) factorizes as: 

n , M , fp(\p\) 8p(z) Po n n 

r \p\ l \+z z \p\ 

The denominator |p| 2 (l+z 2 ) explicitly accounts for dimensions. Without loss of generality we can 
choose gp(0) = 1. The data for gp(z) = (1 + z 2 )Dp(\p\,p )Dp(\p\,0)~ l are shown in Fig. [jjfor 
L = 24. Their leading behaviour can be well described by a power law (1 + z 2 ) a - For j8 > 2.3 the 
functions gp vary consistently with L, violating multiplicative renormalizability. However as L — > oo 
all values of a are compatible with 1 within one or two a i.e. Dp(\p\,po) might eventually be po 
independent. Eq. ( |2.1| ) has deep consequences on the "naive" calculation of D(\p\) oe Y, PQ D(\p\,po) 
as in ^ 0, ^ . Consider different lattice cut-offs for space and time ^ = £ > 1 and define p = a s \p\. 
Neglecting subleading terms, for large L we approximate the sum over po by an integral yielding: 

Dp{p )^ / 2 ^-Dp(p , Po )= H . . / — (l+r) = i /(—«); 
p y_A 2tt p |p| Jo Jt |p| p 

where B(z,a,b) is the incomplete beta function. In the lattice Hamiltonian limit, corresponding 

to § -> oo [JTJ, / becomes |/?| independent 1 , l(^f,a) -> 2 ^r(i-") - Then D j3(I^D " and 
multiplicative renormalizability relies solely on (|/?|). In a standard lattice formulation, however, 
t, = 1 and the extra \p\ dependence S( 4^1, 5, — 05 + j) cannot be avoided. Fig. |2| shows, in spite 



of approximations, nearly perfect agreement between Eq. ( |2.2| ) and the slopes observed in the UV 
for the naive definition of D(\p\) as in 

l I can be analytically continued if a > a — j^N. 
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Figure 1: Data for gn (z) vs 1 + z 2 in log-log scale, L — 24. For sake of readability not all j3 are shown. 



The above discussion makes clear that the static propagator should be defined as Dp(\p\) 
^yfS 1 . To extract it at available L and j8 we fit gp (z) and define: 



1+z 2 ~ 
fp(\p\) = Dp{\p\,po)—-r =■■ Dp(\p\,p ) 



(2.3) 



which is now independent of po, up to noise. To improve the signal we average over po, yielding: 

(2.4) 



M\P\) ■= jZDp(\p\,Po) , Dp(\p\) := 



Fig. H shows the resulting Dp(\p\), which is multiplicatively renormahzable. Fitting a power law 
in the IR, \p\ q M~ l ~ q , and a power law plus logarithmic corrections in the UV, M r ~ x log \p\ \~ s 
gives q = 0.99(1), r = 1.002(3), s = 0.002(2) and M = 0.88(1) GeV, with £ 2 /d.o.f. all in the range 
3.2-3.3, in agreement with the UV and IR analysis in [|], [14], |l^]. We thus constrain q = r = 1, 
s = and fit the whole result through the Gribov formula: 



D{\P\) 



2,/LpI 2 + ^ 



(2.5) 



We find just as good agreement (# 2 /d.o.f. = 3.3) again with M = 0.88(l)GeV ~ 1\[c . 
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Figure 2: Comparison between UV deviations from the Gribov formula in the MC data [j7j] for the naive 
static propagator Y,p Q Dp(\p\,po) (•) and our prediction for the leading term B(j£p, j, — a + j) ( — ). 

3. Conclusions 

We have shown that on the lattice the static transverse gluon propagator is multiplicatively 
renormalizable, IR-UV symmetric and can be well described by Gribov 's formula over the whole 
momentum range. Its infrared and ultraviolet behaviours are in good agreement with the results 
obtained in the variational approach to continuum Yang-Mills theory in Coulomb gauge [Q, O]. 
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Figure 3: The gluon propagator: MC data (•), (a few) data for \p\ — > M 2 \p\~ l (OX the fit to Gribov's 
formula ( — ) and the result of the Hamiltonian approach [|l4j] ( — ). 
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